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ABSTRACT 



Theory and simulations suggest that magnetic fields from radio jets and lobes 
r \ ■ powered by their central super massive black holes can be an important source 

ph ■ of magnetic fields in the galaxy clusters. This is paper II in a series of studies 

P-r where we present self-consistent high-resolution adaptive mesh refinement cos- 

mological magnetohydrodynamic (MHD) simulations that simultaneously follow 
+r£ ■ the formation of a galaxy cluster and evolution of magnetic fields ejected by an 

active galactic nucleus (AGN). We studied 12 different galaxy clusters with virial 
masses ranging from 1 x 10 14 to 2 x 10 15 M Q . In this work we examine the effects 
of the mass and merger history on the final magnetic properties. We find that the 
evolution of magnetic fields is qualitatively similar to those of previous studies. 
In most clusters, the injected magnetic fields can be transported throughout the 
cluster and be further amplified by the intra-cluster medium (ICM) turbulence 
during the cluster formation process with hierarchical mergers, while the amplifi- 
cation history and the magnetic field distribution depend on the cluster formation 
and magnetism history. This can be very different for different clusters. The total 
magnetic energies in these clusters are between 4 x 10 57 and 10 61 erg, which is 
mainly decided by the cluster mass, scaling approximately with the square of the 
total mass. Dynamically older relaxed clusters usually have more magnetic fields 
in their ICM. The dynamically very young clusters may be magnetized weakly 
since there is not enough time for magnetic fields to be amplified. 

Subject headings: galaxies: active - - galaxies: clusters: general - - methods: 
numerical — MHD — magnetic fields 
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Introduction 



The i ncreasing detections of radio emission from galaxy clusters, calle d radio halos and 
relics (see lCarilli Sz Taylorll2002t Ferrari et al.ll2008t iGiovannini et al.ll2009l ) suggest that the 
ICM is permeated with magnetic fields. Radio halos are generally diffuse and extended over 
> 1 Mpc, covering the whole clusters, while radi o relics, which are often o bserved at the edges 
at clusters, can extend as long as 2 Mpc (e.g. Ivan Weeren et al.l 120101 ) . By assuming that 
the total energy in relativistic electrons is comparable to the magnetic energy, the magnetic 
fields in the cluster halos can be estimated at 0.1 — 1.0 /iG and the total magnetic energy 



can be as high as 10 61 erg (lFerettilll999l ). 



Magnetic fields in galaxy clusters are also extensively studied through Faraday rota- 
tion measurements (FRM). Distribution of FRM, combined with the ICM density mea- 
surement, o ften yields cluster mag netic fields of a few to ten /iG level, mostly in the clus- 
ter centers ( jCarilli fc Taylor! 120021 ). More interestingly, FRM was investigated to suggest 
th at cluster magnetic fields may have a Kolmogorov-like turbulent spectrum in core regions 



rtVogt fc Enfilinl 120051). with an energy spectrum p eak at several kpc. Other studies (for ex- 
ample iTaylor fc Perleyl Il993t lEilek &: Owenl |2002| ) have suggested that the coherence scales 
of magnetic fields can range from a few k pc to a few hundred kpc, implying large amounts 



of magnetic ene rgy and fluxes in the ICM ([Colgate fc Li 



200 0J). Recently stud y of magnetic 



fields by FRM ( iGovoni et al.ll2006t iGuidetti et al.ll2008t Bonafede et al.ll2010h are extended 
to the outer part of clusters by observations of more radio galaxies behind or embedded 
in these clusters. It is expected that the Extened Very Large Array (EVLA) will provide 
unprecedented new observations on the magnetic fields in the ICM. 

It is unlikely that magnetic fields have been dynamically important during the cluster 
formation. But, it is suggested that the strength and geometry of magnetic fields in clusters 
may play a crucial role in cluster formation through other processes, such as heat transport, 
which conse quently aff ect the applicability of clusters as sensitive probes for cosmological 
parameters (jVoitll2005l ). In addition, since magnetic fields are closely related to synchrotron 
emission and FRM, the distribution of magnetic fields is important to the understanding the 
radio observations of the ICM. 

Magnetic field evolution is highly nonlinear during cluster formation and difficult to be 
studied analytically. Cosmological MHD simulations are being used to study the proper- 
ties of magnetic fields in the ICM. Such simulations can be very useful in interpreting the 
magnetic field information f rom the observ ation results. Since the origin of cluster magnetic 
fields is still being debated ( IWidrowl 120021 ). various initial magnetic fields are used in MHD 
cluster formation simulations. Simulations were done with initial magnetic fields either from 



some random or uniform fields at high redshifts (IDolag et al.ll2002l ; iDubois &: Teyssierl 12008 



-3- 



Dubois et al.l 120091 ) or from the outflows of normal galaxies (jDonnert et al.l 120091). All these 



simul ations have found that fields can be further amplified by cluster merger (IRoettiger et al. 



19991 ) and turbulence in addition to collapse, and their findings roughly match results from 
observations. On the other hand, very small se ed fields from some first princ iple mechanism, 



like Biermann b attery effect, are also studi ed (IKulsrud et a 



19971 : Kul 12009 ) in galaxy clus- 



ter simulations. IKulsrud et all ( 119971 ) and iRyu et al.l ( 120081 ) suggested that very week seed 
fields can be amplified by dynamo processes in clusters, though current simulations have not 
been able to model such processes self-consistently. 

Magnetic fields in large scale radio jets and lobes from AGNs serve as one very intriguing 
source of cluster magnetic fields because observa t ions show that they can carry large amounts 



of magnetic energy and fl ux (jBurbidgel Il959l ; iKronberg et al.l l200ll ; ICroston et al.l 12005 



McNamara fc Nulsenl 120071 ). The magnetization of the ICM and the wider inter-galactic 



medium (1GM) by AGNs has been suggested o n the energetic grounds ( iColgate fc Lil 12000 



Furlanetto &; Loebll200ll ; IKronberg et al.ll200ll ). without details of the physical processes of 
magnetic field t ransportation and amplification. Through cosmological MHD simulations, 
Xu et al.l (120091 ) showed that magnetic fields injected from a single AGN into a local region 
can be sufficient to ma gnetize the whole cluster to the micro Gauss level by the ope ration 
of sma ll-scale dynamo (jBrandenburg fc Subramanianl 120051 ; ISubramanian et al.l 120061 ) . Re- 
cently, ISur et al.l ( j2010l ) used MHD simulations to show that same process can also operate 
during the formation of first stars. Small-scale d ynamo may be an im portant process to 
generate magnetic fields in various cosmic objects (jSchleicher et al.ll2010l ). 



In the first paper (jXu et al.ll2010l ) of this study, we studied the magnetic field evolution 
in a single massive, relaxed galaxy cluster with different AGN injected energies and injection 
redshifts. That paper found that, as long as the magnetic fields are injected before the active 
merger period during the cluster formation history, the AGN magnetic fields can be spread 
throughout the whole cluster and get substantial amplification. The final magnetic fields are 
weakly dependent on the amount of initial magnetic fields. But the behavior of magnetic 
fields in clusters of different masses and at various dynamical states, young (unrelaxed) or 
old (relaxed) clusters, is still not well studied. 

In this paper, we perform a series of high resolution adaptive mesh refinement (AMR) 
MHD simulations of 12 galaxy clusters of different masses with initial magnetic fields injected 
by an AGN. This allows us to investigate the robustness of magnetizing the ICM using the 
AGN magnetic fields and address additional questions that are not examined with a single 
cluster. We explore the properties of magnetic field distribution and their evolution in the 
ICM of those clusters in different dynamical states during their formation histories. The 
organization of this paper is as follows. In Section [2J we provide the details of the simulation 
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setup, including the cluster formation and the magnetic fields injection. We then summarize 
the main results in Section |3j We present the detailed spatial distribution of magnetic 
fields, evolution of the magnetic energy and the radial profiles of magnetic fields strength. 
We also present and discuss the properties of synthetic Faraday rotation measurement of our 
simulated clusters. In Section 4, we present a summary of the main findings and conclusions. 



Basic Model and Simulations 



We have performed self-consistent adaptive mesh refinement (AMR) cosmological MHD 
galaxy cluster formation simulations of a set of clusters with initial magnetic fields injected by 
an AGN, using cosmological MHD code ENZO +MHD JCollins et al.lboioh . The simulation 
set up is the same a s that in Ku et al.l (J2010I ). The initial conditions are generated from 



an 



Eisenstein &: Hul (119991 ) power spectrum. The simulations use a ACDM model with 
parameters h = 0.73, Q m = 0.27, f2 fe = 0.044, Q^ = 0.73, and cr s = 0.77 . These parameters 



are close to the values from recent WMAP 3 results ( iSpergel et al.ll2007l ). Initial conditions 
from different random seeds are used to generate different clusters. The simulation volume of 
each run is 256 h~ 1 Mjpc on a side, and it uses a 128 3 root grid and 2 level nested static grids 
in the Lagrangian region where the cluster forms. This gives an effective root grid resolution 
of 512 3 cells (~ 0.69 Mpc) and dark matter particles of mass resolution of 1.07 x 10 10 M o . 
During simulations, 8 levels of refinement are allowed beyond the root grid, for a maximum 
spatial resolution of 7.8125 h~ l kpc. 

The simulations are evolved from redshift z = 30 to z = using an adiabatic equation 
of state, with gamma=5/3. The Simulations do not include other physics, such as radiative 
cooling or star formation feedback, as they are not important to the majority of the dynamics 
of cluster formation. Here, we have simulations of 12 galaxy clusters with m ass ranging from 
9.9 x 10 13 to 1.9 x 10 15 M at z — 0. In previous studies ( IXu et al.ll2010l ). we have shown 
that the magnetic field evolution is neither sensitive to the injection redshift between z=3 
and z=l nor to the injected AGN magnetic energy, so in this study we have only one injected 
AGN magnetic energy and a single injection time. We "turn on" the AGN with magnetic 
field injection at redshift z = 3 in the most massive halo (or the second most massive one in 
some runs). The initial magnetic field s are injec t ed into the I CM locally assuming that they 
are fro m an AGN (se e description in Ku et al.l (120081 120091 )). basing on a magnetic tower 
model ( )Li et al.ll2006l ). The injected magnetic energy is ~ 6 x 10 59 erg for all runs. Based 
on the dynamical states of these simulated clusters, these 12 clusters fall into 2 groups at 
z=0, each has 6 clusters. We consider the clusters, which have more than half of their final 
masses by z=0.5, as relaxed clusters, and the clusters, which gain more than half of their 
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final masses from z=0.5 to z=0, as unrelaxed clusters. Group one is the relaxed clusters. 
They are labeled as Rl to R6, in the order of their final masses. Group two is the unrelaxed 
clusters. They are labeled as Ul to U6, also in the order of their final masses. For two 
of the unrelaxed clusters (Ul and U2) we have performed two simulations, with magnetic 
energy injected into different progenitor halos. They are labeled as Ula, Ulb, U2a and U2b, 
respectively. The reason that we have injections in different locations is because the choice 
of initial magnetized halo turns out to play an important role in the final magnetism of these 
two clusters. So we have a total of 14 simulations. The properties of the injection halos and 
the final clusters are summarized in Tabled], while the injection magnetic energies are listed 
in Table H 



The AMR settings are the same as in Ku et al.l (J2010I ). The AMR is applied only 
in a region of (~ 50 Mpc) 3 where the galaxy cluster forms. During the course of cluster 
formation before the magnetic fields are injected, the refinement is controlled by baryon 
and dark matter overdensity. After magnetic field injection, all the regions where magnetic 
field strength is higher than 5 x 10~ 8 G are refined to the highest level, with 7.81 h -1 kpc 
resolu tion. The data analysis in this paper was performed using yto for Enzo (JTurk et al. 



201lh . 



There are limitations with our current simulations, which should be keep in mind when 
interpreting the results from these simulations. The highest resolution of ~10 kpc, which 
already made these simulations very big and very computationally expensive, is close to the 
characteristic scale of the ICM magnetic fields obtained from observations. So our simula- 
tions miss some small scale features of the magnetic fields. Higher resolution simulations are 
underway. In addition, the MHD treatment without any kinetic effects in our simulations are 
not enough to completely understand the magnetic properties in the ICM. In some situations, 
the difference introduced by kinetic effects may be dramatic. For example, when effects of 
anisotropic pressure is include d in a weakly-collisional plasma, growth rate of magneti c field 



strength may be much higher (ISchekochihin fc Cowleyll2006l ; ISchekochihin et al.ll2008[ ). Un- 



fortunately, the microphysical processes are still impossible to be simulated self-consistently 
in cosmological simulations with current computational capability. 



1 http://yt.enzotools.org 
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Table 1: Properties of simulated galaxy clusters. 





Cluster Properties 


(z=0) 


Cluster Properties 


(z=3) 


Name 


r 200 (Mpc) 


M vir (M Q ) 


M gas (M Q ) 


r 20 o(Mpc) 


M vir (M Q ) 


M gas {M Q ) 


Rl 


2.161 


1.252el5 


1.863el4 


0.194 


1.577el3 


2.296el2 


R2 


1.909 


8.633el4 


1.215el4 


0.171 


l.lllel3 


1.495el2 


R3 


1.590 


4.985el4 


7.116el3 


0.160 


9.132el2 


1.223el2 


R4 


1.364 


3.149el4 


5.000el3 


0.191 


1.548el3 


2.109el2 


R5 


1.148 


1.877el4 


2.864el3 


0.190 


1.536el3 


2.160el2 


R6 


0.927 


9.897el3 


1.491el3 


0.171 


1.113el3 


1.541el2 


Ula 


2.498 


1.934el5 


2.661el4 


0.191 


1.552el3 


2.278el2 


Ulb 


2.468 


1.866el5 


2.733el4 


0.143 


6.540el2 


9.612ell 


U2a 


1.743 


6.572el4 


8.697el3 


0.089 


1.578el2 


1.940ell 


U2b 


1.718 


6.292el4 


7.963el3 


0.142 


6.349el2 


9.033ell 


U3 


1.709 


6.198el4 


8.105el3 


0.116 


3.460el2 


4.269ell 


U4 


1.675 


5.829el4 


7.665el3 


0.126 


4.449el2 


5.954ell 


U5 


1.626 


5.338el4 


7.370el3 


0.133 


5.318el2 


7.007ell 


U6 


1.447 


3.763el4 


5.312el3 


0.146 


6.908el2 


9.070ell 
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3. Results 

3.1. Formation of the Galaxy Clusters and Evolution of Magnetic Field 

Distribution 

We first briefly show the formation histories of our simulated clusters. The hierarchy 
formation histories of our simulated clusters are presented by plotting the projected gas 
densities at various redshifts from z=3 to z=0. Relaxed clusters are plotted in Figure [TJ 
and the unrelaxed clusters are shown in Figure |5J The clusters usually undergo numerous 
mergers at redshift between 3 and 0.5. For redshift less than ~ 0.5, the clusters start to relax 
unless some rare big mergers happen. The dynamically old and relaxed clusters in group 1 
have experienced major mergers several Gyr before the end of the simulation, and have been 
dynamically relaxed for some time before z=0. In contrast, the younger, unrelaxed clusters 
of Group 2 have experienced major mergers closer to the end of the simulations. As we 
have mentioned at previous section, the clusters in group 1 have more than half of their final 
masses at z=0.5, while the clusters in group 2 more than double their masses from z=0.5 to 
z=0 by big mergers. From the shapes of the projected density, some clusters of group 2 look 
relaxed , e.g. cluster Ul, though they have just finished a very big merger and their ICM 
motions are still quite active. The cluster U2 is undergoing active mergers after z=0.25 and 
is the dynamically youngest cluster in our simulations. This makes the evolution of magnetic 
fields in U2a and U2b very different from other runs. 

Figure [3] and H] show the projected magnetic energy densities. Magnetic fields mostly 
follow the motion of the halos in which magnetic fields are injected and may move a long 
distance during the cluster formation. No matter where the magnetic fields are injected, 
the final magnetic fields tend to distribute around the cluster centers at z=0. These plots 
clearly show that mergers, especially major mergers, spread the magnetic fields throughout 
the clusters. So the distribution of magnetic fields are effected by when these big mergers 
happen during the cluster formation history. There is a clear trend that the sizes of the 
magnetized areas and their associated magnetic energy densities are proportional to the 
sizes of their host halos. This is no surprise since larger clusters, which are formed by more 
mergers and have higher turbulence level to diffuse, amplify and maintain their magnetic 
fields. 

Magnetic fields have already been spread to a large area of the whole clusters when 
simulations stop at z = in all runs, except for runs U2a and U2b. But the distribution 
of magnetic fields can be quite different from that of the ICM gas if major cluster mergers 
have occurred recently. In such case, the magnetic field strength may peak away from the 
cluster centers. Magnetic fields in larger clusters, which are less affected by late time small 



mergers, are distributed more regularly. 

In most runs, the halo with the initial magnetic injection grows fast and becomes the 
major one progenitor of the clusters during their formation histories, so their magnetic field 
evolution and distribution are similar within their groups. But that is not true in three 
clusters, R2, Ul and U2. For cluster R2, though it is well formed before z = 0.5 and has 
relaxed at final output, its magnetized halo does't merge with the major one until z « 1. 
So its magnetic fields can't spread to a large volume and get amplified in the early time. 
This makes its magnetic fields weak and only locally distributed as a large relaxed cluster. 
For the two unrelaxed clusters Ul and U2, they don't have a single major halo during their 
cluster formation until very low redshift. Initial magnetic fields are then injected into two 
different halos to see the different magnetic filed evolution. The cluster Ul is finally formed 
by the merging of two similar size sub-clusters at z ~ 0.2. In runs Ula and Ulb, we initially 
magnetize one of these two sub-clusters, respectively. Their late major merger doesn't have 
enough time to spread their magnetic fields to the other parts of the cluster. The cluster U2 
is formed by active mergers of several sub-clusters at very low redshifts. Runs U2a and U2b 
have two of them magnetized initially. In both runs, the lack of big mergers in the early time 
causes the magnetic fields to have much less amplification, and their final magnetic fields 
are very weak and locally distributed at z=0. The magnetic fields in run U2b are so weak 
that they are barely seen in Figure H] 



3.2. Magnetic Energy Evolution 

We present the evolution of the total magnetic energy of all the simulations in Figure 
separated into two groups of different dynamical states. The magnetic energy generally 
decreases rapidly initially, for a few hundred million years, because of the rapid expansion 
of the initial magnetic structure. It typically takes about 1-2 Gyr for the magnetic fields to 
expand in the ICM to a large volume to catch the turbulent motions and to be amplified. 
Since we inject large amount magnetic energy into the system, which may be much more 
than the magnetic energy that some small halos, like in cluster R6, can maintain. In such 
cases, their magnetic energies are keep decreasing until they drop to the levels that their 
magnetic fields can be sustained in their host halos by the turbulence. After that the 
magnetic energy gradually increases due to the increasing size and the ICM turbulent level 
by continuous mergers until saturation occurs (due to lack of mergers) or the simulation 
ends. Although the precise evolution of magnetic fields is quite different for various runs, 
the fastest amplifications of magnetic fields generally happen a while (hundred million years) 
after the active mergers (big or continuous mergers) which assemble the major part of the 
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clusters. 



The final magnetic energies are between 4 x 10 and 1 x 10 bl erg. Th e magnetic 



energy is saturated (stops growing) (for a detailed discussion, see lXu et al.ll2010l ) for several 
Gry for relaxed clusters in group 1. The only exception is cluster R2, whose weak magnetic 
fields, due to lack of early major mergers, are still being amplified long after the relaxation 
of the cluster. Almost all of the clusters in group 2 have major mergers at low z, so their 
magnetic energies keep increasing until simulations end, except for run U2b. For this run, 
its magnetic fields, which are in a very small halo before the final assembly, occupy too small 
a volume to get amplified effectively when the simulation ends. 

Not surprisingly, the total magnetic energy is higher for the larger clusters, which have 
deeper potential well, and higher thermal and kinetic energies. We plot the total magnetic 
energy versus the virial mass of the clusters in Figure [6] at z — 0. The total magnetic energy 
is scaled as Ml iral for most of the runs, except runs U2a and U2b. The total magnetic energy 
and the virial mass can be a single proportionality parameter a as : 



E 



M 



6.09 x 10 60 er# 



o 



lVl v i r 

1O 15 M ' 



where a is within a range between 0.131 to 1.445 for all runs except runs U2a and U2b. 
The energy 6.09 x 10 60 is chosen to make a = 1 for run Rl. We summarize a, magnetic 
energies and their ratio to kinetic energies inside R200 of all simulations in Table HJ The 
coefficient a is mostly decided by the cluster formation history. They are generally higher 
for the relaxed clusters in group 1 than for the dynamically younger clusters in group 2. 
The magnetic energies of runs U2a and U2b of the youngest cluster are far away from the 
E^j ~ M 2 scaling. This is because their magnetic fields do not get amplified until the late 
active mergers, and their current magnetic energies are limited by their small halos before 
the major active mergers. 

Though the processes of cluster formation and magnetic field evolution are highly non- 
linear and complex, this scaling relation between magnetic energy and virial mass can be 
qualitatively understood. Using an overly- simplified isother mal model, the thermal and ki- 
netic energies of galaxy clusters are proportional to M vir ( Bryan fc Normanl Il998l ) with a 
radial density profil e as p oc r~ 2 and a flat ve locity dispersion. These relations were numer- 
ically confirmed by IVazza et al.l (120061 . 1201 ll ) using both SPH and AMR codes. Since the 
clusters are not exactly isothermal, especially in the outer part of the clusters, these scalings 
should be modified somewhat. The profiles of the densit y and the velocity dispersion deviate 
from the isothermal distribution as p oc r -3 and a oc r (INavarro et al.lll995t ISunyaev et al. 



20031 ) at large radii. In such a case, we find that the total kinetic energy is approximately 
proportional to the square of the virial masses in our simulations. As the magnetic energies 
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in clusters in similar dynamical state are proportional to the kinetic energies with small 
scattering related to formation and magnetization histories (since magnetic power spectra 
of all simulations have a same shape and their levels are simply determined by the levels of 
the kinetic ones (see Section. I3.4p ). the magnetic energy is then proportional to the square 
of the cluster mass. 

Since the magnetic energy in a halo is proportional to the square of it mass, the magnetic 
fields in a massive cluster are much larger than the magnetic fields residing in its progenitor 
halos. So the major contributing factor to the final cluster magnetic fields is from the dynamo 
process after majority of the cluster has formed. It also suggests that additional magnetic 
fields from more AGNs or regular galaxies in smaller halos may have only small impact on 
the final magnetic fields in a cluster. 

As listed in Table [2J the magnetic energy in a cluster is only a small fraction of its 
kinetic energy. For relaxed clusters, their magnetic energies are about 1% of their kinetic 
energies inside their virial radii. These ratios drop to about 0.1% for the dynamical younger 
clusters, and even smaller for the recently formed cluster U2. The local magnetic energy 
is much smaller than the kinetic energy in most of the cluster regions as well. To show 
this, we plot the volume histogram of kinetic /3 {gk/^m) of all clusters at z=0 in Figure 
[7J For relaxed clusters, kinetic /3 peaks at between 50 and 100, while peaks of kinetic (3 
are bigger than 100 for unrelaxed clusters. Even for a relaxed cluster, only a small fraction 
of its volume has magnetic energy bigger than kinetic energy. So, the magnetic fields are 
dynamically unimportant in all these simulated clusters. Several factors may contribute to 
why the magnetic energy is only a small fraction of the corresponding kinetic energy. First, 
it is possible that simulation resolution is not high enough, as the corresponding numerical 
diffusion is big, so the turbulence decays before it has an opportunity to amplify the magnetic 
fields. Some of our higher resolution test runs show that there are somewhat larger amount 
of more magnetic fields, but the total magnetic energy is still much smaller than the kinetic 
energy. Second, it could be that there is not enough time during cluster formation for 
the magnetic fields to grow in large scales, of which the eddy turnover time is ~ Gyr. 
Third, the driving from mergers is not constant and not homogeneous, and becomes weaker 
when the cluster gets bigger. Fourth, it is not known what the expected level of magnetic 
field saturation is for compressibl e turbulence. Simula tions of supersonic and superalfvenic 



compressible MHD turbulence in iKritsuk et al.l (120091 ) found that the saturated magnetic 



energy is much smaller than the kinetic energy if the initial fields are weak. A further study 
to understand why the magnetic fields saturate at a level much lower than the kinetic energy 
in the ICM is needed. Fifth, since we only have one seed source for magnetic fields, significant 
magnetic fields only fill a portion of the ICM (from 20% to 80% see later section on magnetic 
fields distribution). So in a case of more sources and a higher filling of magnetic fields in 
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a galaxy cluster, its magnetic energy ratio could potentially be somewhat higher (but still 
much smaller that kinetic energy, at most a few percent). 



3.3. Magnetic Field Distribution over the ICM at z = 

3.3.1. Radial Profiles of Magnetic Fields 

In Figure [SJ we present the spherically rms averaged radial profiles of magnetic fields at 
z=0. The relaxed and unrelaxed clusters are plotted in two different panels, and the radii 
are normalized by their virial radii. It is clear that the profiles of magnetic field strength are 
determined both by the cluster sizes and their dynamical states. 

The magnetic field strength is generally higher for the bigger clusters and/or dynamically 
older clusters. For relaxed cluster, micro Gauss fields are present at the centers of clusters 
bigger that 10 15 M , but fields drop to about 0.5 fiG for the smallest cluster of 10 14 M . The 
magnetic field profiles of relaxed clusters in group 1 are more regular as a function of radius, 
while the magnetic field profiles of dynamically younger clusters in group 2 are obviously 
disturbed by their recent big mergers. For runs U2a and U2b, since the cluster is formed 
so late and no significant amplification of magnetic fields has been operated, their magnetic 
fields are much weaker and reside only in the inner part of the cluster. 



Several other MHD simulations (e.g. iDolag et al.l 120021 ; iDubois et al.l 120091). a s well as 



the study on the RM and the X-ray surface brightness correlation (JDolag et al 



20011 ) . suggest 



that the decline of magnetic fields is correlated with the gas density as n 7 . Analysis of RM 
and X-ray brightness data in lDolag et al.l (|200ll ) yields the 7 = 0.9 for A119, and 7 = 0.5 for 
3C129 with large uncertainty. We fit our magnetic field radial profile with the gas density 
profile and list the best-fitted results in Table [2] . The 7 in our simulations scatters between 



0.3 and 0.81, which is consistent with the results from IDolag et al.l ( 120011 ). Relaxed clusters 
usually have smaller 7, so have flatter magnetic field profiles as they have more time to 
amplify their magnetic fields in the outer part by the ICM turbulence after major mergers. 
However, the reliability of our results may be limited by the way we magnetized clusters, 
we need simulations with many AGNs in one cluster to get more reliable n 7 profiles. In 
addition, since this result is obtained by fitting the averaged radial profiles of the magnetic 
field strength and the gas electron density, it doesn't mean that there is a simple correlation 
between the magnetic field strength and the gas density throughout the cluster. 
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Table 2: Properties of magnetic energy in simulated galaxy clusters 



Name 


E\ f (erg) a 


El , (erg) b 


a c 


771200 / 771200(2 


7 e 


Rl 


6.17e59 


9.55e60 


1.00 


0.00814 


0.44 


R2 


4.85e59 


1.34e60 


0.295 


0.00306 


0.65 


R3 


5.57e59 


9.10e59 


0.601 


0.00240 


0.57 


R4 


5.32e59 


4.19e59 


0.693 


0.0122 


0.3 


R5 


6.37e59 


3.07e59 


1.436 


0.0113 


0.46 


R6 


5.94e59 


8.62e58 


1.445 


0.0113 


0.45 


Ula 


6.32e59 


2.99e60 


0.131 


0.00119 


0.56 


Ulb 


5.83e59 


5.88e60 


0.277 


0.00248 


0.55 


U2a 


3.32e59 


8.90e58 


3.38e-2 


1.77e-4 


0.67 


U2b 


6.01e59 


3.68e57 


1.53e-3 


7.19e-6 


0.81 


U3 


4.74e59 


4.09e59 


0.175 


8.21e-4 


0.43 


U4 


5.35e59 


1.25e60 


0.603 


0.00513 


0.74 


U5 


5.35e59 


6.07e59 


0.350 


0.00224 


0.59 


U6 


5.21e59 


3.29e59 


0.382 


0.00259 


0.53 



a: magnetic energy at the end of injection; 

b: magnetic energy inside R 20 o at z=0; 

c: a in Equation [H proportionality parameter between magnetic energy and cluster mass 

square; 

d: ratio of magnetic to kinetic energy inside R200 a t z=0; 

e: best-fit of 7 in relation between the magnetic field strength and the gas density (B oc n 7 ); 
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3.3.2. Magnetic Field Spatial Distributions 

We plot the volume histograms and complementary cumulative volume ratio histograms 
of magnetic field strength inside the virial radii for all simulations at z=0 in Fig. |9j The 
magnetic field fillings are quite different for the relax and unrelaxed clusters. For dynamically 
older clusters, most volumes inside the virial radii are filled with magnetic fields. Magnetic 
fields stronger than 1CT 8 G typically fill more than 75% of the cluster volume for group 
one except run R2, and less than 50% for group two except run U5. For run R2 of group 
one, since the magnetized halo merges with the bigger halo after z=l, its magnetic field 
volume filling is consequently smaller. Once the magnetic fields are well spread throughout 
the clusters and get amplified, magnetic field strength distribution peaks between 0.1 and 
0.3 /iG. Larger and older clusters have much higher peaks as well as longer tails of stronger 
magnetic fields. Runs U2a and U2b have very weak magnetic fields that only fill very small 
volumes. 

We also plot the two-dimensional distribution of the magnetic field strength verses the 
gas density inside their virial radii in Figure [TUl showing how magnetic fields distributed 
over the ICM plasma. There is no obvious correlation between the field strength and the 
gas density. The distribution is similar for clusters of similar sizes. Most of the magnetic 
field strength is between 0.1 to a few micro Gauss, and are mixed with gas over a wide range 
of densities. This casts doubts on cluster magnetic field modeling when simple correlation 
between \B\ and ICM density is assumed. 



3.4. Kinetic and Magnetic Energy Density Power Spectra 

The kinetic and magnetic energy density power spectra of all clusters at z=0 are shown 
in Figure dU which are computed from boxes of 512 3 cells (~(5.5 Mpc) 3 ) in the highest level 
enclosing the clusters. The ICM turbulence is represented by Kolmogorov-like spectra in 
kinetic energy. These kinetic spec tra are also seen from pure hydrodynamics simulations of 



galaxy clusters ( iVazza et al.ll2009l ). All magnetic spectra are in a similar shape and follow the 



k 3 / 2 Kazantsev law in the large scales. The Kazantsev spectrum, which is the prediction o f 



small-scale dynamo theory in the kinematic regime ( see iBrandenburg fc Subramanian 



Federrath et a 



is also found in simulations of a collapsing Bonnor-Ebert sphere in 

Thes e magnetic spectra in our simulations show that the small-scale dynamo (see 



2005) 



(1201 lh. 



Xu et al. 



20101 . for detailed discussions) does operate in all our simulated clusters. For dynamically 
young clusters, some (like cluster U4 ) of their power spectra are disturbed by the recently 
mergers. 
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For the relaxed or dynamically older clusters, the magnetic fields have enough time to 
be amplified, their magnetic energy densities are close to the kinetic energy densities in small 
scales. One possible reason that the magnetic energies in smaller scales are smaller than the 
kinetic energies is because the magnetic fields do not fill all the space of the clusters. For 
the recently merged clusters, the magnetic fields do not have enough time to spread through 
out the newly formed clusters and be completely amplified by the ICM turbulence. So their 
magnetic energies are much smaller than the kinetic ones even in the small scales. 



3.5. Faraday Rotation Measurement from Simulations 



The Faraday Rotation Measurement is a key method for measuring cosmic magnetic 
fields and has provided important information on cluster magnetic fields. Th e RM maps have 
been u sed not only to estimate the strength of magnetic fields (see review by Carilli &: Taylor 



(120021) 1 but also the turbulent structure of the fields (IVogt &: Enfilml 12005 



Enfilin fc Vogt 



20061 1 in the ICM. We compute the RM maps using the magnetic field and free electron 
density distribution from our simulations by integrating 8 Mpc over the clusters centered at 
each cluster's center along the y axis at z=0. We show our synthetic RM maps in Figure 
[T2l T hough the basic featu res of the RM distribution are similar to those from previous 
study ( iXu et al.ll2009l . 120101 ). clusters of different size and/or at different dynamical states 
have quite different RM distribution in their absolute values, cover areas, and structure 
scales. The morphologies of the RM maps reflect the dynamical states of the clusters and 
their merger histories. There is a clear trend that there are more small scale patterns in 
the dynamically old relaxed clusters than in young clusters. This is because they have more 
time for the magnetic field amplification in small scales after major mergers. On the other 
hand, there are more long filaments on the younger, unrelaxed clusters reflecting their recent 
mergers. There are very long filaments (> 1 Mpc) in runs Ula and Ulb associated with 
their recent large head-on mergers. It is interesting to see that the RM maps are quite 
different between these two runs, of which the pre-merger magnetic fields reside in different 
sub-clusters, though the gas dynamics of these two runs are almost identical. This shows 
that the magnetism history of the cluster also plays an important role in determining their 
RM distribution. 

The 2-D azimuthally averaged radial profiles of the absolute values of these FRM (|i?M|) 
are plotted in Figure [T3J while the radial profiles of the standard deviation of the RM are 
shown in Figure HH The \RM\ profiles resemble the magnetic field strength profiles with 



steep slopes, as the ICM gas density decreases with the r adius. These \ RM \ profiles are 
roughly similar to the pattern from observational data in IClarke et al.l fJ200l[ ). The <7rm 
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profiles are also consistent with the recently observational results in iGovoni et al.l ( 120101 ) 



We also plot the area histograms of RM and the complementary cumulative histograms 
of \RM\ inside 500 kpc central circles of the clusters in Figure US Relaxed clusters generally 
have more areas covered by significant |i?M|. For all cases, the distribution of positive 
and negative RM is roughly symmetric. This is because there is no net magnetic fluxes in 
the clusters, for the model of field injection we have used. The RM area histograms are 
similar for clusters of similar sizes and dynamical states, though their RM maps may look 
quite different from their different formation or magnetism history. This suggests that the 
histograms of RM distribution are not sensitive to the cluster merger and magnetism history, 
which can be reflected by their RM filament and patchy structures. 



4. Conclusions 

In this paper, we report an ensemble of simulations of magnetic field evolution in galaxy 
clusters with a wide range of masses between 1 x 10 14 and 2 x 10 15 M in various dynamical 
states. With similar amounts of initial magnetic fields injected from a single AGN at a 
high redshift, all clusters are eventually filled with micro Gauss magnetic fields, except for 
dynamically very young clusters. The power spectra of kinetic energy density show that 
the ICM is in a turbulent state, while the spectra of magnetic energy density show that the 
small-scale dynamo process is being driven by by the ICM turbulence. This result, along with 
the previous study of a single cluster with the magnetic field injections of different amounts 



of magnetic energy and at different redshifts (IXu et al.l 120101 ). suggests that magnetization 



of galaxy clusters by the operation of small-scale dynamo with the seed magnetic fields from 
AGN is very robust, and it produces magnetic fields consistent with observed magnetic fields. 

The magnetic field evolution and distribution are decided both by the masses of the 
clusters and their dynamical formation histories. The total magnetic energy is approximately 
scaled as the square of the virial mass of the cluster, while the dynamically older (relaxed) 
clusters usually have more magnetic fields. This implies that the cluster magnetic fields 
are mostly determined by the dynamo process of the ICM turbulence generated by the 
hierarchical mergers. Additional magnetic fields from more AGNs or smaller cluster forming 
halos will not have a major impact on the final magnetic fields in a massive cluster. The 7 in 
the scaling relation between magnetic field and gas density radial profiles (\B\ oc n 7 ) range 
between 0.3 and 0.81 for our simulated clusters, while relaxed clusters have flatter magnetic 
fields profiles. In addition, the relaxed clusters usually have self-similar magnetic field radial 
distribution, while the field distribution in younger unrelaxed clusters is disturbed by their 
recently mergers. Though our simulated clusters only have initial magnetic fields from a 



-16- 



single local source, most volumes in the simulated clusters are well magnetized. 

We also studied Faraday rotation measurements obtained from the magnetic field and 
gas density distribution in our simulated clusters. They are also determined by both the 
cluster sizes and their dynamical states. The radial profiles of \RM\ resemble the profiles of 
the magnetic field strength, and are consistent with the pattern from the observational data. 
The RM maps reflect the recent cluster mergers, as well as the cluster magnetism history. 
There are very long filaments in the RM maps in the recently merged clusters, while their 
small scale patchy bands reflect the ICM turbulence. Very different distributions of RM of 
the same cluster but with magnetic fields injected in different locations are observed. This 
suggests that RM distribution may be a good probe not only for the cluster formation but 
also its magnetism history. A detailed study on the relation of RM from simulated clusters 
with their ICM turbulence will be presented in a forthcoming paper. 

This work was supported by the LDRD and IGPP programs at LANL and by DOE 
office of science via CMSO. Computations were performed using the institutional computing 
resources at LANL. ENZO_MHD is developed at the Laboratory for Computational Astro- 
physics, UCSD with partial support from NSF AST-0808184 to M.L.N. 
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Fig. 1. — Projected baryon density for relaxed clusters in group 1 at different redshifts. The 
first column is projections at the time of injection, and each image covers a region of 32.0 
Mpc x 32.0 Mpc comoving. The second column shows at z=2 and each image size is also 
32.0 Mpc x 32.0 Mpc comoving. The third, fourth and fifth columns are images at z=l, 
z=0.5, and z=0.25, respectively. The plot sizes are 16.0 Mpc x 16.0 Mpc comoving. The 
final column shows projections at z=0. Each image covers an area of 8.0 Mpc x 8.0 Mpc. 
The first panel in each row is marked with the respective runs. The color range of all plots 

g cm" 



is from 1.0 x 10~ 4 to 1.0 x 10" 2 - — " 2 
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Fig. 2— Same as Figure. (TJ but for unrelaxed clusters in group 2. 
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Fig. 3. — Projected magnetic energy density for relaxed clusters at different redshifts. The 
first column shows at z=2 and each image size is 32.0 Mpc x 32.0 Mpc comoving. The 
second, third and third columns are images at z=l, z=0.5, and z=0.25, respectively. The 
plot sizes are f 6.0 Mpc x f6.0 Mpc comoving. The final column shows projections at z=0. 
Each image covers an area of 8.0 Mpc x 8.0 Mpc. The first panel in each row is marked 
with the respective runs. The color range is from 1.0 x 10 8 to 1.0 x 10 12 erg cm -2 . 
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Fig. 4.— - Same as Figure. but for unrelaxed clusters in group 2. 
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Fig. 5. — Evolution of total magnetic energy of the simulations, the upper panel shows 
relaxed clusters, while the bottom panel shows the unrelaxed clusters . 
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Fig. 6. — Total magnetic energy vs. the virial mass of all simulations at z=0. 
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Fig. 7. — Volume histogram of kinetic plasma j3 (= &k/&m) inside the clusters' virial radii 
at z = 0. Relax clusters are shown in the top panel and the unrelaxed clusters are shown in 
the bottom panel. 
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Fig. 8. — Spherically averaged radial profiles of rms magnetic field strength of different 
clusters. Relax clusters are shown in the upper panel and the unrelaxed clusters are shown 
in the bottom panel. The x-axis is normalized by their virial radii. 
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Fig. 9.— Volume distribution (top) and complementary cumulative volume distribution 
(bottom) of the magnetic field strength inside the virial radius. The left panel shows the 
relaxed clusters, while the right panel shows the unrelaxed clusters. 
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Fig. 10. — Volume weighted two-dimensional distributions of the magnetic field strength vs. 
baryon density of all runs inside the cluster virial radiuses at z = 0. Contour lines are the 
volume of gas at that density and magnetic field at 10 fc cm3, where k = 69.0, 69.2, 69.4, 
. . . 71.2 from outer to inner. The upper two rows are relaxed clusters, and the rest are 
unrelaxed clusters. 
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Fig. 11. — Kinetic and magnetic energy density spectra of relaxed (Top) and unrelaxed 
(Bottom) clusters. Kinetic and magnetic energy densities of same cluster are represented 
by same line. Kinetic energy densities of all clusters show k~3 Kolmogorov-like spectra. 

3 

All magnetic energy densities have a short ka Kazantsev spectrum, while the levels of the 
magnetic energy are depend on the cluster dynamical states and their magnetized histories. 
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Fig. 12. — Faraday rotation measurement of the clusters by integrating through the cluster 
on the y direction. It covers a region of 3 Mpc x 3 Mpc at z = 0. The color range shown is 
from —500 (blue) to 500 (red) rad m -2 . The upper two rows are relaxed clusters, and the 
rest are unrelaxed clusters. 
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Fig. 13. — Azimuthally averaged radial profiles of \RM\ of the RM maps shown in Figure 
T2l The top panel shows the relaxed clusters, while the bottom panel shows that unrelaxed 
clusters. The x-axis is normalized by the clusters' virial radii. 



-33- 



10 c 



„10' 



CO 



DC 

10 1 



10 L 






0.6 0.8 

radius [R v ] 



R1 
R2 
R3 
R4 
R5 
R6 




1.2 



10 c 




0.6 0.8 

radius [R v ] 



U1a 

U1b 

U2a 

U2b 

U3 

U4 

U5 

U6 



1.2 



Fig. 14. — Azimuthally averaged radial profiles of standard deviation of RM of the RM 
maps shown in Figure [T2j The top panel shows the relaxed clusters, while the bottom panel 
shows that unrelaxed clusters. The x-axis is normalized by the clusters' virial radii. 
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Fig. 15.— Histogram of RM (top) and complementary cumulative histogram of \RM\ 
(bottom) of the central circle of 0.5 Mpc radius of the RM maps shown in Figure [T2J The 
left panel shows the relaxed clusters, while the right panel shows the unrelaxed clusters. 



